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PREFACE 


The  author  thanks  Dr.  A.  Bratenahl  for  pointing  out  prior  appearances 


of  the  present  magnetospheric  model  in  the  literature. 
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INTRODUCTION 


A standard  exercise  in  magnetospheric  plasma  physics  [e.  g. , 

Roederer,  1970]  is  to  calculate  the  location  of  the  plasmapause  and 
of  other  equipotencial  trajectories  for  a given  model  of  the  electric 
(E)  and  magnetic  (B)  fields.  Since  it  is  customary  to  assume  that 
E ■ B = 0 and  that  9E/9t  = 0,  one  can  derive  E = - VV  from  the 
scalar  potential  V = V(o,  p),  where  a and  P are  the  Euler  potentials 
f Stern,  1967,  1971]  from  which  one  obtains  the  magnetic  field 
JB  = ( Vo  ) X.  (Vj  P )•  In  other  words,  the  electrostatic  potential 
depends  only  on  coordinates  that  remain  constant  along  magnetic 
field  lines.  Thus,  an  equipotential  trajectory  comprises  the  family 
of  field  lines  visited  by  a particle  drifting  in  the  direction  of  E X B. 

There  have  appeared  several  calculations  of  such  (cold-plasma) 
drift  trajectories,  but  only  the  dipolar  model  of  ^ [e.g.,  Alfven. 

1939;  Roederer,  1970;  Grebowsky,  1970,  1971;  Kivelson  and  Southwood, 

1975;  Cowley  and  Ashour-Abdalla,  1976]  has  thus  far  yielded  analytical 

expressions  for  these  equipotential  trajectories.  More  realistic  models  j 

of  B^  [e.g.,  Taylor  and  Hones,  1965;  Taylor,  1966;  Brice,  1967; 

Kavanagh  et  al. , 1968]  have  yielded  jE  X^.B  trajectories  only  in  the  form 
of  numerical  output.  The  purpose  of  this  letter  is  to  report  on  the 
analytical  calculation  of  cold-plasma  drift  trajectories  in  a certain 
non-dipolar  model  of  B>,  i.  e.  , in  the  field  model  obtained  by  super- 
imposing a uniform  magnetic  field  parallel  to  the  magnetic-dipole  axis  j 

of  the  planet  [e.g.,  Dungey,  1961;  Hill  and  Rasabach,  1975],  The  resulting 


"magnetosphere"  has  a tail  that  is  coaxial  with  the  dipole  but  topo- 
logically equivalent  to  the  infinite  geomagnetic  tail  idealized  in 
various  magnetospheric  models  [e.  g.  , Ness,  1969].  Moreover,  the 
present  model  should  apply  equally  well  to  the  known  magnetospheres 
of  other  planets. 


BASIC  EQUATIONS 

The  magnetic  field  in  the  present  model  is  given  [Hill  and 
Rassbach,  1975] 

B = B_  a V[ ( a/r  )2  cos  6 - (a/b)?  (r/a)  cos  0],  (li 

AMV 

where  Bq  is  the  equatorial  (0=  tt/2)  magnitude  of  the  dipolar  cont-i- 
bution  to  B at  the  planetary  surface  (r  = a)  and  b is  a constant  that 
ultimately  characterizes  the  magnetospheric  diameter.  The  spherical 
coordinates  r,  0,  and  f are  measured  respectively  from  the  center 
of  the  planet,  the  north  pole,  and  the  midnight  meridian.  The  equation 
of  a magnetic-field  line  is  found  to  be 

r = [l  + (1  /2)(r /b)3]  La  sin2  0 (2) 

2 

if  one  requires  that  sin  0 - r/La  as  r — 0.  Since  (1)  is  derivable 
[ Stern,  1967]  from  the  Euler  potentials 

o = - (a2/r)  Bq  [1  + (l/2)(r/b)3]  sin2  0 (3) 
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and  p = a <p,  the  parameter  L ( = - a Bg/o)  thus  serves  as  am 
invariant  field-line  label. 

5|C 

It  follows  from  (2)  that  some  field  lines  are  open  (L.  > L ) and 
some  are  closed  (L<L  ).  The  configuration  is  illustrated  in  Figure  1, 
The  separatrix  (L  = L ) can  be  identified  from  the  property  that  the 
direction  of  B becomes  indeterminate  around  its  equator,  which  is  to 
say  that  B = 0 there.  Thus,  it  follows  from  (1)  that  the  separatrix 

7 MA  AM 

crosses  the  equatorial  plane  (0  = tt/2)  at  r = b,  whereupon  it  follows 
from  (2)  that  L = 2b/3a. 

Open  field  lines  (L  >L  ) admit  the  limit  r— co.  It  follows  from 
(2)  for  open  field  lines  that 


p s r sin  0 - (2b3/La)2 /2  = (3L  '7l)] /2  b 


as  z (s  r cos  0)  approaches  infinity.  Thus,  the  magnetospheric 

diameter  2 p (=  2b3  ' j,  attained  asymptotically  as  | z | — co,  follows 

at  once  from  (4).  The  magnetopause  corresponds  to  the  boundary  (L  - L ) 

between  closed  and  open  field  lines.  The  open  field  lines  in  Figure  1 

correspond  to  the  tail  of  the  magnetosphere. 

The  functional  form  of  the  electrostatic  potential  V(o,  (3  ) = 

V(L,  <p)  remains  to  be  specified.  If  the  planet's  rotation  at  frequency 

fl/2n  is  coaxial  with  its  "magnetosphere"  in  the  present  model,  the 

corresponding  ^-independent  term  in  V(L,  <p  ) is  given  [e.  g. , Schulz, 

2 

1970]  by  - O a B^/cL.  No  other  form  is  acceptable.  The  ^-dependent 
term  representing  magnetospheric  convection  must  be  continuous  at 


* 

L = L . 


It  is  convenient  to  take 


r 


* 

H 


V (L,  ip  ) = Ej  (L/L*)n  La  sin  <f  - (ft  a2  BQ/cL)  (5  a) 

(where  n is  a constant)  for  l < jl  and 

V(L,V)  = (4L*/9L)1/2  Ejb  sin  <p  - (na2BQ/cL)  (5b) 

for  L > L . The  first  term  of  (5a)  corresponds  to  an  electrostatic 
potential  of  the  general  form  postulated  by  Volland  [1975], 

The  first  term  of  (5b)  corresponds  to  an  asymptotically  uniform 
electric  field  across  the  magnetospheric  tail  [Hill  and  Rassbach,  19751. 

This  is  a special  case  of  the  electrostatic  potential  postulated  by 

Volland  [ 1975]  for  L > L*.  However,  Volland  [1975]  assumed  a 

* 

dipolar  magnetic  field  B and  was  thus  required  to  introduce  an  L 

quite  arbitrarily  in  (5).  The  present  discontinuity  in  E^  at  L - L it  i 

an  unavoidable  consequence  of  the  present  magnetic  field  model 
which  was  introduced  by  Dungey  [1961]  and  resurrected  by  Hill  and 
Rassbach  [1975]. 


TAIL  CONVECTION  PATTERN 

$ 1 /z 

Since  (4)  indicates  that  p — • (3L  /L)  b as  | z | — * <d  , the  quantity 

£ 1 f? 

(L  /L)  in  (5b)  should  serve  as  an  expedient  quasi  radial  coordinate  in 
the  tail  of  the  magnetosphere.  The  E^X.B.drift  °*  c°ld  plasma  there 
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r 


generates  a family  of  constant-V  surfaces  that  corresponds  to  a family 

$ 1/2 

of  coaxial  cylinders  if  (L  /L)  is  plotted  as  a function  of  y>  in  polar 
coordinates.  The  cylinders  have  their  common  axis  at 


,1/2 


*2 


(L  /L)  ' = j cEj  L *72  O aB0  = (1/2)  (L  / Lj ) 


(6) 


on  the  dawn-dusk  meridian.  The  axis  has  an  invariant  longitude 
f = (tt/2)  sgn  ( f2a/c)  if  (as  is  assumed  here)  Ej /BQ  > 0.  The  equation  of 
any  such  constant-V  cylinder  is 

[(L  /L)  ^ sin  <?  - (1/2)  (L  / Lj  )2  sgn  (tta/c)]2 

+ (L*/L)  cos2  ? = (1/4)  (L  /Lj)^  - (cL*V/  SI  a^B^),  (7) 

and  the  right-hand  side  of  (7)  represents  the  square  of  the  radius. 

. * -1/2 

If  Lj  / L >2  , then  it  follows  from  (6)  that  the  common 

axis  of  the  drift  cylinders  must  lie  within  the  magnetotail.  Suffic  iem’v 
small  drift  cylinders  in  this  case  will  lie  entirely  on  open  field  lines, 
and  will  thus  preclude  the  adiabatic  access  of  solar-wind  plasma  to 
their  interiors.  Drift  cylinders  that  bear  labels  V such  that 


- (V/E1L"a)sgn(na/c)  < (Lj/L*)2-! 


(8) 


have  this  property,  and  the  boundary  of  the  adiabatically  inaccessible 
cylinder  just  grazes  the  magnetopause  at  <p  = (w/2)  sgn  (Sla/c  ) . 

The  foregoing  results  are  illustraced  by  quantitative  example  in 

>!< 

Figure  2 for  three  distinct  values  of  Lj /L  . Magnetotail  convection 
of  the  correct  sense,  i.  e.  , from  cos  <p  < 0 toward  cos  » > 0,  is 
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W 


I~  Mercury) 


(Earth) 


('•Jupiter) 


Figure  2.  Asymptotic  (z  = ±<jo)  cross  sections  of  the  tail  of  present 
magnetosphe ric  model  showing  cold-plasma  flow  patterns  ( fl  > 0)  tor 

jf  * 1/2 

three  selected  values  of  L,  /L  . Radial  coordinate:  (L  /L) 


assured  by  the  assumption  (see  above)  that  Ej/Bq  > 0.  Cases  for 
which  Lj/L  < 2 provide  adiabatic  access  of  cold  solar-wind 

plasma  to  all  the  tail  field  lines,  since  the  corresponding  drift  shells 
are  partial  cylinders  (see  Figure  2b).  The  concentric  drift  cylinders 

a(c 

degenerate  to  parallel  planes  in  the  limit  Lj/L  = 0 (see  Figure  2a). 

* -1/2 

Cases  for  which  Lj/L  > 2 lead  to  the  adiabatic  exclusion  of  cold 

solar-wind  plasma  from  a cylindrical  region  of  the  tail  adjacent  to  the 
dawn  (dusk)  meridian  for  > 0 (ft  < 0),  as  illustrated  in  Figure  2c. 


OMEGAPAUSE  LOCATION 

Conventional  magnetosphe ric  models  follow  Nishida  [l  966]  in 
identifying  the  plasmapause  as  the  boundary  between  closed  and  open 
equipotentials  of  (5a).  Such  an  identification,  however,  entails  no 
estimate  of  the  actual  plasma -density  profile  (which  depends  also  on 
the  rate  of  plasma  transport  along  13,  for  example).  Thus,  Siscoe  and 
Chen  [l  975]  have  suggested  "omegapause"  as  an  appropriate  name  for 
the  boundary  between  closed  and  open  equipotentials  of  (5a).  The  value 
of  V on  the  omegapause  in  a given  field  model  is  usually  found  by  locating 
a stagnation  point  in  the  equatorial  convection  pattern,  i.e.,  a closed 
field  line  on  which  W s 0.  This  corresponds  to  a field  line  on  which  the 
direction  of  E X EMs  arnl:,iguous . It  follows  from  (5a)  that  the  coordinates 
of  this  field  line  must  be  f = - (tt/Z)  sgn  ( fta/c)  and 

L = flaB0/cE1L5!<2(n  + 1)  j 1 + 

= [ L1Z/(n  + l)L**]1/(n  + 2)L*.  (9) 
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where  Lj  = ( fiaB^/cEj  | *^2  as  in  (6).  If  /L*  < (n  + l)*^2,  one 
obtains  the  express  ion  [cf.  Voiland,  1 97 5 ] 


(L/L*)n+*  sgn(  fia/c)  sin  <p 


(Lj/L*)2  (L*/L)  - {n  + 2)[L2 /(n+ I)L*2](n+ I,/(n  + 2) 


(10) 


for  the  trajectory  of  the  omegapause.  This  form  allows  one  to  plot 

$ 

without  difficulty  as  a function  of  L/L  , as  in  Figure  3b.  The  omega- 
pause  bears  the  label 

2 , , . ,v  . *2  T (n+  l)/(n+2) 


V = 


(n  + 2)  E L a [L,  / (n  + 1)  L ] 


sgn(S2a/c),  (11) 


which  is  the  value  of  the  electrostatic  potential  to  which  it  corresponds. 

jjj  1/2 

The  case  L^/L  > (n  + 1)  yields  no  closed  field  line  on  which 
E = 0.  In  this  case  the  largest  closed  equipotential  surface  must  be 
the  one  that  just  grazes  L = L at  <p  = - (tt/2)  sgn  ( fla/c).  This 
grazing  trajectory  bears  the  label 


= - Ejl/a  [l  + (1^/L*)2]  sgn  ( Sla/c) 


(12) 


and  generates  an  omegapause  profile  given  by 
(L/L*)n+  1 sgn  ( f2a/c)  sin  ip 

= <VL*)2  [(L*/L)  - 1]  - 1 (13) 

A representative  profile  of  this  type  is  plotted  in  Figure  3c.  Equi- 
potential trajectories  exterior  to  (13)  necessarily  extend  onto  open  field 
lines  (L  > L*)  and  thus  presumably  lose  their  plasma  of  ionospheric 
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Figure  3.  Equatorial  (z  = 0)  cross  sections  of  the  present  magneto- 
spheric  model,  showing  locations  of  the  omegapause  for  selected 
values  of  Lj/L  and  of  n.  Omegapause  is  degenerate  in  (a),  the 
case  of  a nonrotating  planet.  Radial  coordinate:  l/l". 


origin.  Brice  and  Ioannidis  [l970]  had  qualitatively  sketched  (for  Jupiter) 
an  omegapause  that  just  grazed  the  magnetopause  at  dawn.  The  present 
equations  illustrate  such  a phenomenon  in  closed  analytical  form. 

The  mathematical  condition  for  having  an  omegasphere  that 
grazes  the  magnetopause  is  L^/L  > (n  + 1 ) . The  condition  for  the 

magnetosphe ric  tail  to  contain  closed  convection  trajectories  (in  the  sense 
of  Figure  2c)  is  Lj/L*  > 2-1^2.  The  similarity  of  these  two  conditions 
makes  it  unlikely  that  one  of  them  will  be  satisfied  without  the  other 


also  being  satisfied.  Only  for  (n  + 1) 


r (n  + 1)1/2  > L , / L > Z 


there  be  a eusped  omegapause  on  closed  field  lines  (Figure  3c, 
solid  curve)  in  the  same  planetary  magnetosphere  with  open  (tail) 
field  lines  that  are  adiabatically  isolated  (as  in  Figure  2c)  from  the 
magnetospheric  surface  (L  = L ).  This  transition  (1  j/L  ~ 0.7-1. 5) 
between  convection-dominated  ( L j / I-  « 1 ) and  rotation -dominated 

(Lj/L  » 1)  magnetospheres  spans  only  a factor  ~2  in  I j/L 
if  one  takes  n = 1 as  recommended  by  Volland  [1975]  . The  case  of 
a nonrotating  planet  (Figures  2a  and  3a)  corresponds  to  Lj/L  =0 
and  represents  an  extreme  case  of  the  former  topology.  In  this  case 
the  omegasphere  shrinks  to  zero  radius  and  disappears  beneath  the 
planetary  surface. 

Application  of  the  foregoing  remarks  to  the  plasmapause  is 
contingent  on  the  identification  of  it  with  the  "last"  closed  electrostatic 
equipotential,  i.e.,  with  the  omegapause,  as  in  Nishida  [1966].  Such 
an  identification  is  contingent  on  the  rate  of  plasma  transport  along  J3, 
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as  Siscoe  and  Chen  [1975]  have  noted.  It  is  also  contingent 

on  the  lack  of  plasma  transport  across  equipotentials  of  (5a).  Lemaire 

[1974,  1975]  has  strongly  challenged  this  latter  assumption  on  which  the 

identification  of  the  plasmapause  with  the  omegapause  is  contingent,  and 

has  argued  that  interchange  motions  driven  by  the  effective  gravitational 

field  exert  a controlling  influence  on  the  shape  of  the  plasmasphere.  It 

would  be  fair  to  say  (with  respect  to  this  issue)  that  the  jury  is  still  out. 

$ 

Thus,  it  should  be  emphasized  that  the  criteria  on  L^/L  derived  here 

relate  to  the  shape  of  the  omegasphere  and  to  the  magnetotail  convection 

* 1/2 

pattern.  Whether  the  inequality  between  Lj/L  and  (n  + 1)  ultimately 
determines  the  shape  of  the  plasmasphere  is  a question  that  remains 
unanswered  in  the  light  of  reservations  noted  above. 


INTERPLANETARY  SCALING 


Insertion  of  the  usual  terrestrial  parameters  (a  = 6371.2  km,  Bq  = 

0.308  gauss,  E^  = 0.  6 V /km,  ft  = 2ir  day  *)  in  (6),  with  appropriate 

I 1 1 /2 

attention  to  physical  units,  yields  L^  = ftaBQ/cEj|  = 4.877  irre- 

j{e 

spective  of  n.  The  choice  of  L = 8.  547  leads  to  open 

field  lines  at  all  invariant  latitudes  above  70°,  to  an  asymptotic  tail 
3 -5 

field  B^a/b)  = 14.76  X 10  gauss,  and  to  an  asymptotic  diameter 
2p  = 44.37a  for  each  lobe  of  the  geomagnetic  tail.  All  these  parameters 
are  quite  reasonable  and  correspond  to  a value  of  0.571  for  the  ratio 
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Lj  /L  . Thus,  the  earth's  magnetosphere  corresponds  to  Figures  2b  and  3b. 
It  would  be  of  great  interest  in  the  present  context  to  generate 

jjc 

estimates  for  Lj /L  in  the  magnetospheres  of  Jupiter  and  Mercury. 

If  calculated  from  the  model  of  Mendillo  and  Papagiannis  [1971],  the 

parameter  |E^|  scales  as  (a^lB^R^)  ^ ^ for  a constant-velocity 

solar  wind  f Schulz,  1975],  where  R is  the  distance  from  the  sun.  The 

value  of  L * should  scale  as  ( | Bq  | R)*  ^ , i.  e.  , as  the  distance  from  the 

center  of  the  planet  to  the  magnetopause  (measured  in  planetary  radii). 

* 231 /4 

Thus,  the  ratio  L^/L,  should  scale  as  ( | B^  ) a ) . Insertion  of 

5jC 

accepted  planetary  parameters  in  this  scaling  law  yields  L^/L  ~ 10 
for  Jupiter  and  Lj/L*~  10-^  for  Mercury.  The  latter  result  suggests 
that  Mercury  (with  L.  ~ 2)  has  no  omegasphere,  since  it  implies  that 
Lj^  1.  Mercury's  magnetosphere  must  resemble  that  illustrated  in 
Figures  2a  and  3a. 

The  estimate  of  L^/L  ~ 10  suggests  that  the  magnetosphere  of 
Jupiter  is  a more  extreme  example  of  the  case  illustrated  in  Figure  2c. 

For  Jupiter  the  offset  cylinder  of  drift  shells  adiabatically  isolated  from 
the  solar  wind  represents  about  99%  of  the  tail  volume.  The  accessible 
crescent  contains  only  about  1%  of  the  tail  volume.  Conversely,  the 
oval  omegasphere  covers  about  98%  of  the  equatorial  plane  of  Jupiter's 
magnetosphere  (excluding  the  tail).  The  remaining  crescent,  adiabatically 
accessible  to  the  solar  wind,  covers  only  about  2%  of  the  equatorial  area. 

The  foregoing  results  suggest  that  solar-wind  plasma  has  no  access 
to  most  of  Jupiter's  magnetotail  except  by  diffusion  across  the  adiabatic 
(£.£ trajectories.  This  mode  of  access  should  be  rather  inefficient. 
Thus,  the  present  calculation  has  shown  how  the  access  of  solar -wind 
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plasma  to  most  of  Jupiter's  magnetotail  is  at  least  inhibited  by  planetary 
rotation.  One  might  conclude,  other  things  (such  as  1 B^  a3  > being 
equal,  that  a nonrotating  Jupiter  would  have  had  a much  denser  tail 
plasma.  Alternatively,  one  might  look  to  planets  such  as  Jupiter  (with 

jjc 

large  /L  ) for  examples  of  tail  plasma  derived  largely  from  the 
planetary  ionosphere  [cf.  Axford,  1970]  rather  than  from  the  solar  wind. 

However,  the  present  results  also  have  major  implications  for  the 
earth's  magnetosphere.  Although  it  is  found  (see  Figure  2b)  that  solar- 
wind  plasma  has  adiabatic  access  to  field  lines  across  the  entire 
geomagnetic  tail,  the  adiabatic  trajectories  themselves  show  a 
pronounced  curvature  in  consequence  of  the  earth' s rotation.  It  might 
be  possible  to  verify  this  curvature  experimentally  by  tracking  the  plasma 
component  of  a barium  cloud  released  there. 

The  present  calculation  applies  to  an  arbitrary  planet,  with  Mercury, 
Earth, and  Jupiter  taken  as  specific  examples.  The  calculation  illustrates 
the  importance  of  superimposing  the  convection  and  corotation  electric 
fields  throughout  the  magnetosphere  (in  the  tail  as  well  as  on  closed  field 
lines).  It  also  emphasizes  the  operational  effectiveness  of  thinking  in 
terms  of  E and  B fields,  rather  than  in  terms  of  "moving"  magnetic- 
field  lines.  The  superposition  of  convection  and  corotation  electric  fields 
is  undoubtedly  implicit  in  the  earlier  work  of  Taylor  and  Hones  [19f>5], 
but  its  consequences  there  are  more  difficult  to  identify  because  the 
output  was  purely  numerical.  The  present  calculation,  based  on  the 
magnetospheric  model  of  Hill  and  Rassbach  f 1975],  enables  all  the  physical 


boundaries  and  limiting  trajectories  to  be  identified  in  terms  of  precise 


analytical  expressions.  Of  course,  the  present  results  would  be 
distorted  somewhat  in  their  application  to  a more  realistic  magneto- 
sphere, but  the  basic  topology  of  plasma  drift  in  the  tail  of  a rotating 
planetary  magnetosphere  should  correspond  well  to  that  found  here. 

i! 
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